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ABSTRACT 

Near-infrared adaptive optics imaging with NICI and NaCO reveal what ap- 
pears to be a three- winged or lobed pattern, the "butterfly nebula" , outlined by 
bright Br7 and H2 emission and light scattered by dust. In contrast, the [Fe II] 
emission does not follow the outline of the wings, but shows an extended bipolar 
distribution which is tracing the Little Homunculus ejected in r] Car's second 
or lesser eruption in the 1890's. Proper motions measured from the combined 
NICI and NaCO images together with radial velocities show that the knots and 
filaments that define the bright rims of the butterfly were ejected at two different 
epochs corresponding approximately to the great eruption and the second erup- 
tion. Most of the material is spatially distributed 10° to 20° above and below 
the equatorial plane apparently behind the Little Homunculus and the larger SE 
lobe. The equatorial debris either has a wide opening angle or the clumps were 
ejected at different latitudes relative to the plane. The butterfly is not a coherent 
physical structure or equatorial torus but spatially separate clumps and filaments 
ejected at different times, and now 2000 to 4000 AU from the star. 

Subject headings: circumstellar matter - ISM:individual(Homunculus Nebula) - 
ISM:jets and outflows - stars:individual(?7 Carinae) - stars:mass loss - stars:winds, 
outflows 
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Introduction 



Eta Carinae's dusty, mottled, bipolar Homunculus nebula and its ragged equatorial 
"skirt" of debris provide a visual record of its two eruptions in the IQth century. Numerous 
studies to measure the proper motions and expansion of the ejecta have shown that the 
bipolar lobes were ejected during rj Car's "great eruption" in the 1840's (Gaviola 1950, 
Ringuelet 1958, Currie et al. 1996, Smith & Gehrz 1998, Morse et al. 2001). Proper 
motion and radial velocity measurements of the equatorial ejecta revealed material very 



likely expelled in two separate eruptions (jPavidson et al 



Currie fc Dowling 



1999 



Morse et al. 



1997 



2001 



Smith fc Gehrz 



1998 



20011), in the 1840's and in the lesser known second 
eruption in the 1890's (see Humphreys, Davidson and Smith 1999). A second, smaller 
bipolar outflow, now called the "Little Homunculus", was discovered by Ishibashi et al. 
(2003) based on the Doppler-shifted morphology of the integral-sign shaped Fe H and [Fe H] 
emission lines in HST/STIS spectra. It is embedded in the Homunculus and extends ~ 2" 



in either direction along t he maior axis 



with the 1890's eruption (jlshibashi et al 



ts kinematics suggest an election date associated 



2003 



Ishibashi 



2005 



Smith 



20051). The Little 



Homunculus cannot be distinguished in visual images of rj Car due to bright nebulosity 
near the star and to extinction in the Homunculus, but its structure can be traced by the 
velocities of the [Fe II] emission lines (Ishibashi et al. 2003; Smith 2002, 2005; Teodoro et 
al. 2008). 

Infrared imaging and spectroscopy allow us to penetrate the dusty outer shell of the 
Homunculus to reveal a complex structure near the central star. Mid-infrared imaging, 
from 5/im to 24/im, revealed an elongated feature several arcseconds across centered 
on the star and interpreted as a circumstellar ring (Rigaut & Gehring 1995; Smith et 



al. 1998; Poloms ki et al. 1999: Smith k Gehrz 2000). 



imaging to date, (ISmith et al. 



1995 



2002 



Chesneau et al. 



he hi ghest resolution mid-IR 



20051 ). show what looks like two 
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interlocking rings or loops with multiple bright knots. Near-IR images (1-2 fim) obtained 
with HST/NICMOS were discussed by Smith and Gehrz (1998) who identified a loop-like 
structure which they also attributed to an equatorial torus. A near-IR (1 - 4fim) study by 
Chesneau et al (2005) using the VLT adaptive optics system, NaCCo resolved this feature 
into a complex dusty structure, the "butterfiy" nebula, with three fans or wings outlined 
by emitting warm dust and refiected hydrogen emission, around the relatively dark wings 
which are apparently devoid of dust. The wings of the butterfiy appear to be emanating 
from the region of the central star. Chesneau et al. suggested that this was a dusty, bipolar 
nebula distinct from the little Homunculus. Alternatively, Smith (2006) has argued for a 
distorted equatorial torus that has been swept back by the stellar wind. Kinematic data are 
needed to determine the geometric orientation of this complex structure with respect to the 
axis of symmetry of the Homunculus and the equatorial ejecta and its ejection epoch. 

We obtained narrow-band images with the Near-Infrared Coronagraphic Imager (NICI) 
and adaptive-optics on the Gemini-South telescope just after r] Car's 2009 spectroscopic 
event. In the next section we describe the observing procedure and data reduction. 
The narrow-band Br7, H2, and [Fell] images are discussed in §3. The Br7 and H2 
images are essentially identical. Both trace the outlines of the butterfiy nebula, while 
the continuum-subtracted [Fe II] image reveals a bipolar distribution due to the little 
Homunculus. In §4 we compare our images with the NaCO images obtained several years 
earlier and measure the projected motions. The expansion, orientation, and ejection epoch 
for the resolved features in the butterfiy nebula are discussed in the last section. 



^The Nasmyth Adaptive Optics System and Near-Infrared Imager and Spectrograph on 
UT4 of the VLT under programs 074.D-0140, 076.D-0586, 078.D-0562. 
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2. Observations and Data Reduction 

We used the Near Infrared Coronagraphic Imager (NICI) (Ftaclas, Martm & Toomey, 
2003, Chun et al. 2008) to image the Homunculus of 1] Car on 2009 February, 11 UT, a few 
weeks after the spectroscopic minimum. NICI features a dual channel camera, primarily 
designed for Spectral Differential Imaging (SDI) for planet searches (Marois et al., 2005, 
Biller et al., 2004), but the simultaneity of the two channel imaging also provides optimal 
on/offline imaging. This simultaneity ensures that the observing conditions (transmission, 
seeing) for the two channels are matched. Because of its extended nebulosity, the SDI 
technique was not appropriate for rj Car. 

Pairs of exposures were recorded simultaneously in the two channels with different 
filters. The Cassegrain rotator was turned off during the observations. The filters 
were chosen to map the spatial extent of previously identified emission features in t] 
Car's near-infrared sp ectrum and to provide of f-hand continuum images for subtraction. 



Near-infrared spectra (jSmith fc Davidson 



200ll ) were used to determine which filters would 



yield Oj^-band continuum images free of strong emission features. The filters and their 
characteristics are listed in Table 1. The exposure in each filter is a coaddition of thirty 
0'.'38 s subexposures. During our observations the seeing varied between 0'.'6 and I'.'O. The 
effective resolution was 0'.'056 and 0'.'060 for the H and i^"-band filters, respectively, close 
to the corresponding diffraction limits of 0'.'044 and 0'.'057. NICI has an O'.'OIS pixel scale 
providing a Nyquist sampling of the H and K-hand diffraction limits. We used the 95% 
pupil mask with the 0'.'32 semi-transparent focal plane mask to block the central star. The 
mask slightly degrades the resolution, but minimizes the intensity of the speckle pattern. 
With the rotator turned off, the field slowly rotated (6-8°) through the observing sequence. 
The field rotation was used for spatial dithering and for bad pixel removal. This approach 
is more efficient than a classical dithering pattern as it avoids the overheads involved in the 
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reacquisition behind the coronagraphic mask 

Data reduction steps followed a relatively classical approach. The images were first 
dark-subtracted and flat-fielded using a median-combined fiat that was produced using fiat 
frames taken without the coronagraphic mask. The images were then registered using the 
unsaturated core of the central star, visible through the semi-transparent coronagraphic 
mask. A "bad pixel" mask was produced from the fiat field images, and bad pixels were 
flagged in all images. Images were then derotated to a common position angle using the 
known parallactic angle at the time of observation. The images were median combined to 
produce the final image. This procedure was repeated for the on and off band filters and 
the resulting images were subtracted to produce a line image. Figure 1 shows the resulting 
continuum and continuum-subtracted line images for each of the three filter combinations. 
The combined color continuum and line-continuum images are shown in Figures 2a and 2b. 

To better characterize small scale structures in the images, a set of final images was 
produced by subtracting an 11 x 11 median filter from all of the images before the final 
median combination. These high-passed images better highlight the filamentary structures 
in the inner Homunculus and are used later to measure their motions. 

The near-infrared NaCO images used in this paper were obtained between 2002 and 
2006 by Chesneau and his collaborators. A new detector (InSb Aladdin 3 array) was 
installed in May 2004, replacing the Aladdin 2 detector that had been used since the camera 
was first available (Chesneau et al. 2005). The NB_374, NB_405 narrow band filter were 
used to achieve the best optimization and contrast for the dusty neighborhood of rj Car. 
The field of view was 28" x 28" and the scale was 27.1 mas per pixel. The reading mode was 
the so-called 'Uncorr' mode in which the array is reset and then read once. The minimum 



^Central wavelengths 3.74yum and 4.05/im and 0.02/xm wide. 
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effective integration time used for these observations was 0.1750 s. The Auto Jitter mode 
was used so that at each exposure, the telescope moves according to a random pattern in 
a 10" box. The images were reduced with a self-developed IDL routine that processes the 
individual frames as described in Chesneau et al. 2005. 

3. Description of the Images — Penetrating the Homunculus Nebula 

Our K-hand continuum image and continuum subtracted narrow band Br7 and H2 
images in Figure 1 and the color continuum image in Figure 2a clearly show the three 
winged dark pattern named the butterfly nebula by Chesneau et al. (2005). The dark 
wings in both near and mid-infrared images are attributed to the lack of emitting dust at 
these locations. The continuum-subtracted narrow-band Br7 and H2 images (Fig.l) show 
essentially identical patterns of emission due to scattered light by dust. In Figure 3, we 
show our NICI i^-band continuum image together with the NICMOS narrow-band 2.15//m 
image, at lower resolution, from Smith and Gehrz (2000) and the central region of the 
butterfly from the NaCO and NICI K-hand images. Interestingly, comparison of the two 
high resolution i^T-band images show no apparent change in the morphology of the butterfly 
or of the warm infrared emitting regions. 

In contrast with the Br7 and H2 images, the [Fe II] emission shows a very different 
distribution in Figure 1 and in the combined line color image in Figure 2a. Although the 
corresponding continuum image shows a weak butterfly pattern, the [Fe II] emission does 
not follow the outline of the wings, but instead shows a much broader bipolar distribution 
that is most easily seen in Figure 2b. The spatial extent of the [Fe II] emission is readily 
traced in the nearer SE lobe out to ~ 2" from the star. We also note that a mottled 
structure is detectable in the SE lobe [Fe II] emission probably due to reflection in the 
overlying dusty Homunculus. A relatively bright emission region just above the mask is 
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probably from the Weigelt knots where the [Fe II] emission is known to be strong (Davidson 
et al. 1995, 1997, Smith 2002). Directly to the NW of the star we are very likely detecting 
a contribution to the [Fe II] emission from the "fan" . Despite its appearance in the visual 
images of the Homunculus, this feature is not part of the equatorial spray, but instead a 
region where we see light from the star's north polar region directly reflected off the dust. 
Although fainter in the NW lobe, the [Fe II] emission can also be traced to about 2" along 
the major axis. 

As we noted in the Introduction, the physical extent and shape of the little Homunculus 
i s traced spectroscopi cally by its Fe II and [Fe II] emission. Figure 8 in the discovery paper 



(llshibashi et al. 



20031 ) illustrates the shape and spatial extent of the bipolar [Fe II] emission 



from the Litt 



Homunculus ([Smith 



e Homuncu lus relative to the Homunculus lobes. A velocity map of the 



20061 ) in the 1.64/im [Fe II] line, the same line imaged with NICI, and 
the IFU map by Teodoro et al. (2008) likewise show the extent of the Little Homunculus to 
2" along the major axis. The spectroscopic map also exihibits the same asymmetry in its 
brightness between the SE and NW lobes as mentioned above. We therefore conclude that 
the [Fe II] emission in the NICI image is tracing the Little Homunculus. 

One of our original goals for high-resolution images near the time of rj Car's 2009 
spectroscopic event was to look for possible changes in its dusty environment that may be 
due to changing illumination of the inner ejecta and excitation of the Weigelt knots as the 
secondary star passed through the massive primary's wind. The Weigelt knots are small 
blobs or condensations of gas and d ust 700 AU (0'!3) to the no rthwest of the star first 



identified by speckle interferometry (jWeigelt fc Ebersberger 



19861 ). The high excitation 



emission lines seen in ground-based s pectra originate in t he knots presumably due to UV 



radiation from the hot secondary star flMehner et al. 



20101). During the spectroscopic events 



these emission lines weaken and disappear. 
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As already noted there is no obvious variation in the i^-band images obtained nearly 
seven years apart. Due to the focal plane mask, the region of the Weigelt knots, where 
we might expect to observe some changes, is unfortunately not visible in the NICI images 
although the extremely IR-bright central star can be seen through the mask in the 



ir 



interferometr y (IWeigelt et al. 



spectroscopy ( jPavidson et al. 



19951 



19971 : 



1 

iFalcke et al. 


1996 


) and later in 


HST i 


Borland et al 


2004: 


Smith et al. 


2004) 



moving very slowly, compared with the expansion of the Homunculus, at less than 50 km 
s~^ (0'.'005 yr^^) in the equatorial plane. Fortunately no mask was used for the five 3.74/im 
(Pf7) and 4.05/im (Bra) NaCO images obtained from 2002 to 2006 and which bracket the 
2003.5 event. The structures seen in the near-infrared images closely correspond with the 
knots seen in the visible or UV although they are not spatially coincident. The visible 
structures, dominated by scattering, trace the walls of the dense clumps of dust, while the 
infrared structures are identified with the emission from hot dust, probably the external 
layers of the clumps. The infrared condensations corresponding to knots C and D, are 
therefore called C' and D', see Chesnaeu et al. (2005). 

We carefully measured the photocenter of the two most prominent knots, C' and D' in 
the 3.74/im images at position angles 300° and 352°, respectively, alt hough the extractio n 



20061) 



of the centers was hampered by the increasing brightness of the star (iMartin et al.l 
and Strehl ratio variations. The cuts through the clumps showing the variation of intensity 
vs distance at different times are in Figure 4. The normalized infrared flux appears to be 
decreasing from 2002 to 2005, although we emphasize that accurate photometry is difficult 
with the adaptive optics system. The decrease does not appear to be correlated with the 
2003.5 event or orbital phase. One possibility for the apparent decrease in the IR flux is 
enhanced dust destruction in response to the increased stellar flux, although this will have 
to be verified. 
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The interpretation of the complex structure of the butterfly, with its large asymmetrical 
regions outlined by IR bright emission embedded in the Homunculus, is not straightforward. 
Chesneau et al. (2005) initially suggested that the emission rims of the butterfly wings are 
representative of an outflow along the polar axis projected onto the plane of the sky. Their 
8.7//m image shows bright rims and "dark" regions that correspond approximately with 
the wings in the NaCO images which they ascribe to minimal dust emission. Smith (2006) 
however demonstrated that the H2 spectroscopic emission that outlines the thin walls of the 
Homunculus does not penetrate to the star but appears to terminate at the IR bright edges 
of the dark structure seen in his 8.8//m image. He therefore proposed that this structure is 
actually a disrupted equatorial torus in which the current equatorial stellar wind is clearing 
the equatorial plane of dust. Its complicated shape was then attributed to interaction with 
dusty clumps of different densities. 

To further understand this peculiar structure, kinematic data is needed to measure or 
set limits on its expansion, age, and orientation within the Homunculus. In the next section 
we combine our NICI images with the comparable NaCO images obtained six to seven 
years earlier, to measure the projected motions of the knots and filaments that outline the 
butterfly. 

4. Measurement of the Projected Motions 

We registered all of the NaCO images (K^, 4.05 and 3.74//m) to the NICI K-band 
continuum image by first rotating the images to a common orientation and performing a 
first-order scaling based on the known relative plate scales of NICI and NaCO. The initial, 
relative positioning was done using the unsaturated core of the central star's psf. The 
precise scale and rotation were then fine-tuned iterativcly by maximizing the correlation 
of the high-pass images in an annular region between 2'.'3 and 3'.'5 from the star. This 
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region contains structures that are visible in all of the NaCO and NICI images, but is 
sufficiently far from the butterfly to avoid the inner regions where we want to look for 
anomalous motions. Figure 5 shows the scale vs the date of the NaCO and NICI images. 
Using the scaling factors and the known pixel scales of NACO and NICI, we then derived 
the expansion rate and time since the eruption as a check on our procedure. Using data in 
the range 2'.'3 to 3'.'5 in all of the NaCO and NICI images, we find that the Homunculus 
lobes were ejected in 1850.0 ± 4.0. Using different image combinations, NaCO 3.74/im with 
NICI, NaCO 4.05/im with NICI, and the NaCO data only, we obtain similar results within 
la of the above value. This result is consistent with rj Car's great eruption (1838-1858) and 
previous determinations of the expansion of the lobes cited in §1. 

Using the registered high-pass NICI and the 2002.9 NaCO Ks image, we cross-correlated 
small regions. Since a sharp-edged sample region is very unsuitable for this operation, we 
employed instead a local Gaussian mask with an radius of 0.2". We measured motions 
relative to the co-moving grid that expands with the lobes, i.e., a value of zero corresponds 
to ejection date 1850. The cross-correlation was done in both spatial directions {x,y) with 
shifts ranging from —0.5" to +0.5". The x,y width of each computed cell in the correlation 
datacube was 0.107". We use "cell" rather than "pixel" to avoid confusion with the original 
image pixels which were much smaller, see §2. Because of the Gaussian mask mentioned 
above, the cells are fuzzy and adjoining values are not mutually independent. 

Figure 6 shows the resulting relative co-moving proper motioij^l vector field. At each 
location the correlation amplitude, —1 to +1, indicates the validity of the result; a motion is 
considered valid if the maximum in correlation is greater then 0.4. As expected, the highest 
correlations coincide with strong spatial fine structures, and the periphery of the butterfly 



^Hereafter, we use motion or proper motion with the understanding that these are relative 
motions. 
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has many correlations above 0.8. We experimented with different pairings with the NaCO 
Ks image and obtained similar results for the high significance correlations. The ~ 6 yr 
interval between the two images and the ±0.5" explored by the correlations pose an upper 
limit of ~ 0.01"/yr per axis on the proper motions. All of the high-correlation motions 
are safely below this limit; e.g., Figure 7 shows one elongated structure with a significant 
relative motion of ~ 0.008"/yr. The lower limit on what is measureable depends on how 
much structure there is in a given region; for many regions the la limit is 0.1 mas/yr. 
We used Monte Carlo simulations to estimate rms uncertainties, assuming that the error 
budget is dominated by statistical noise in the shallower NaCO images, estimated from the 
pixel-to-pixel standard deviation using a robust sigma measurement. For every correlation 
cell we adopt the rms error in 100 random trials as the local uncertainty in each component 
of the x,y motioiu. These statistical uncertainties are typically 0.1 to 0.2 mas/yr. 

There are numerous knots, blobs and filamentary features in the vector map. Not 
all of them define the bright rims of the butterfly. For many of these we derived specific 
motions referring to the physical structures rather than cells in the x, y grid. To determine 
the motions of the more recognizable structures we used a circular virtual aperture with 
parabolic weight function 1 — (r/i?)^. (Again we must avoid sharp edges.) For each clear 
physical feature in the map, we first used the virtual aperture to measure a local x, y 
centroid based on brightness. Then, centering the virtual aperture at that position, we 
refer to this weighted local set of cells as an aperture sample or just "aperture." For each 
aperture sample we then determined the weighted mean of the motions for the included 
cells. Experimenting with aperture sizes, we found that the results changed little for 
R ~ 0.16"; therefore we adopted a diameter 2R = 0.321" = 3 cell-widths. The dispersions in 



^Caveat: Obviously this does not include systematic errors, which cannot be quantified 
with the available information. Recall, also, that adjoining cells are not independent. 
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the mean motions for the cells in the apertures are larger than the weighted quadratic mean 
of the cells' individual error estimates. This may be due to systematic effects. Consequently, 
we adopted the standard error of the mean motions for the error in each aperture. We 
retained only aperture samples whose weighted- average correlation amplitudes were greater 
than 0.8. 

The results for the 67 high significance apertures including the angular distance of each 
aperture from the star, position angle, and weighted mean relative motion in R.A. and Dec 
in the co-moving reference frame with their standard errors are in the Appendix. Since the 
motions have been measured in the rest-frame of the expanding Homunculus, we have also 
calculated the intrinsic or true motions in the fixed x, y grid, listed in the last two columns 
of the Appendix. The distribution of these apertures on an image of the butterfly is shown 
in Figure 8 with the vectors for the true or intrinsic projected motions. The corresponding 
velocities in the plane of the sky, at a distance of 2.3 kpc for rj Car (Davidson et al 2001, 
Smith 2002), are given in Table 2 with the epochs when the sub-structures, or knots or 
clumps, were ejected assuming uniform motion radially outward from the star. A histogram 
of the ejection epochs from the projected motions is shown in Figure 9. The names or 
designations of the larger regions used in Table 2 and shown in Figure 10 are adopted from 
Chesneau et al.(2005) and Smith et al. (2003). 

Smith (2006) obtained long-sht spectroscopy of the H2 emission hne at 2.1/xm across 
the Homunculus nebula, and three of the slits, SWl, NEl, and "star", overlap the bright 
rims of the butterfly at five different positions shown in Figure 10. To determine the total 
space motion of the knots or sub-structures in these apertures, their spatial orientation, 
and therefore that of the bright rims of the butterfly, we used the corresponding radial 
(Doppler) velocities, measured relative to the Doppler velocity of rj Car, combined with the 
projected velocity from the proper motions at these five positions. The radial velocities 
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for the apertures within the shts, their resulting total velocity relative to the star, the 
orientation angle relative to the plane of the sky, 9, linear distance from the star, and the 
corresponding ejection epoch are all included in Table 2 for the relevant apertures. 

Throughout this paper we assume that the ejecta move radially outward from the 
central star, and that any apparent transverse motions are measurement errors. Angular 
momentum considerations make this a safe assumption; given the observed total speeds, 
off-center ejection in rj Car's binary orbit would lead to transverse motions less than 0.1 
mas yr^^ at r ~ 3000 AU. 



5. Discussion — Expansion, Age, and Distribution of the Inner Ejecta 

It is perhaps not surprising that the ejection epochs cluster around two dates, circa rj 
Car's great eruption (1838 - 1858) and near 1900 corresponding to its second eruption (1887 
- 1895); see Figure 9. Like many previous authors we find a date after 1895. However, 
acceleration at the 10% level is physically possible and may explain this discrepancy as 
noted by Davidson et al. (1997) (see also Smith et al. 2004) |3 

The orientation angle {9) relative to the plane of the sky has a wide range from 
positive to negative. Most of the sub-regions in Table 2, however, are projected away from 
us at angles from +20° to +60°. The equatorial plane (see the schematic diagram in 
Figure 11) is tilted away from us, behind the SE lobe, and toward us, projected in front 
of the NW lobe, at about 41°. The clumps of knots moving away from us in the SE lobe 

^Note also that in principle the emission center of a condensation could migrate outward 
faster than its center of mass, because the spatial distribution of excitation depends on the 
local UV flux and gas density. This effect is small, however, for any brightness peak that 
represents a real density maximum. 
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thus appear to be distributed above and below the equatorial plane. Either the equatorial 
debris has a fairly wide opening angle of ~ 20° or the material was ejected at somewhat 
different latitudes relative to the plane. Numerous knots or apertures within the larger 
regions have comparable ^'s suggesting that the different projection angles are significant. 
There may also be some correlation between the ejection epoch and the orientation angle, 
although our information for 9 is far from complete. For example, all of the knots with 9 
near 20° (S. Clump and S. Arc) and those with 9 near 50° (SE Clump) have ejection dates 
corresponding to the second eruption, while those consistent with the great eruption have 
orientation angles closer to 40°. Four of the sub-regions, apertures 3, 65, 66, and 67 in 
the "SW Region" and the "NE region", are oriented towards us and are very likely due to 
material projected in front of the NW lobe ejected at the time of the great eruption. 

The total space motions or velocities measured for several of the bright clumps or 
knots range from about 60 to 100 km s~^, much slower that the polar expansion of 
either the Homunculus at 600 km s~^ (see references in the Introduction )lj or the little 
Homunculus at 300 km s^^ (Ishibashi et al. 2003, see also Smith 2005). However, material 
associated with the eq uatorial plane has been found to be expanding at a range of velocities 



( jPavidson et al. 



20011 ) including relatively slow moving ejecta near the star. Davidson et 
al. (2001) also concluded that that the equatorial debris consisted of material from the two 
eruption events. Zethson et al. (1999) reported slow moving gas projected in front of the 
NW lobe with radial velocities relative to the star of —40 to —140 km s~^, and presumably 
in the equatorial plane. The Weigelt knots, mentioned earlier, near the equatorial plane, 
are also slow moving, with radial velocities of ^ —50 km s^^ and an estimated ejection 
epoch of 1900. 



^There have been numerous measurements of the expansion of the Homunculus. See 
review by Davidson and Humphreys (1997) and references therein for the earlier work. 
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Our results for the bright clumps and knots that define the apparent rim of the 
butterfly are thus consistent with results from previous work. They are slow moving, from 
two different eruptions, and near the equatorial plane, but projected away from the plane 
of the sky, 2000 to 4000 AU from the star. In Figure 12 we show the physical location of 
the separate clumps and filamentary features relative to the equatorial plane on a schematic 
diagram of the inner region of the Homunculus. The latitude or angle relative to the 
equatorial plane, (, is also included in Table 2. Given the orientation and distances of the 
sub-regions from the star, most of them are on the far side of the Little Homunculus and are 
also very likely on the far side of the SE lobe. Figure 8 in Ishibashi et al. (2003) llustrates 
the complexity of this region, and also shows a cloud of debris within 1" of the star beyond 
the apparent boundaries of the Little Homunculus and bracketing the equatorial plane, 
similar to the spatial distribution of the sub-regions or clumps in the butterfiy. 

Smith (2006) suggested that the butterfiy is part of a "disrupted torus," related to a 
circular region where the two Homunculus lobes meet near the configuration's midplane. 
This idea can probably be modified to be fairly consistent with the ejecta directions and 
two ejection epochs described above. On the other hand, our results are at least equally 
consistent with other possibilities. Modern investigations (including the work here) reveal 
complex asymmetry, not circular structure, in and around the equatorial plane where the 
lobes would meet in a simple picture. This is not very surprising, since there is no strong 
theoretical reason to expect either circular or latitudinal symmetry there. The central binary 
system is generally acknowledged to have a very eccentric orbit which obviously breaks the 
circular symmetry for outfiows near periastron; in principle the star's equatorial region is 
vulnerable to ins tabilities that may c ause local and chaotic low-speed ejections (e.g., see 



remarks in 55 of ( Zethson et al. 



19991 )). Furthermore, details of the large-scale Homunculus 



lobe shapes an d orientations show that axial symmetry was only an approximation for the 



great eruption (IMorse et al 



1998 



Davidson et al 



200ll ). and large-scale radial "streamers" 
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or "jets" well beyond the lobes have long been known in the images for some time (e.g., 
(Duschl et al. 1995, Meaburn et al. 1987, 1993, 1996). Given these diverse clues, we should 
not expect to represent the past or present equatorial region of the Homunculus in terms 
of two well-defined, intersecting, axially symmetric geometrical lobe-figures. Instead, the 
lobes appear to become ill-defined near the midplane. A large circular configuration around 
7] Car is worth considering, but there is no strong observational or theoretical evidence that 
such a structure exists or existed. Note that the famous "equatorial skirt" - located outside 
the inner region discussed here - consists of radial structures with no discernible circular 
symmetry. 

In the absence of quantitative dynamical models, perhaps the simplest hypothesis is 
that either the unstable primary star can can eject matter asymmetrically, and/or ejections 
can occur near periastron when the companion star's gravity and wind perturb the outflows 
asymmetrically. These phenomena would explain our observations (qualitatively, at least) 
wihout any concentration of pre-ejected material near the midplane. Moreover, in either 
case, local flows can be erratic and asymmetric relative to the midplane because they are 
very sensitive to unstable local conditions. These possibilities were implicit in some of the 
references cited above, but they are extremely difficult to model physically. 

In summary, the appearance of the "butterfly nebula" as a single, coherent structure 
with multiple lobes or wings or as an equatorial torus, was thus due to projection effects 
and to low resolution in the earlier studies. Instead the knots or clumps of knots are 
separate, ejected primarily at two different epochs, and are apparently oriented within ±10° 
to ±25° of the equatorial plane, ~ 2000 to 4000 AU from the star. They may be filamentary 
streamers or jet-like features associated with the equatorial debris on the far side of the SE 
lobe. 

We gratefully acknowledge the support of the Gemini staff and the Director's 
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Pfi. Bendjoya, G. Guerri and A. CoUioud for fruitful discussions and testing of analysis 
methods of the NaCO data. When this work was begun, Artigau was supported by the 
Gemini Observatory, which is operated by the Association of Universities for Research in 
Astronomy, Inc., on behalf of the international Gemini partnership of Argentina, Australia, 
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Table 1. NICI Filter Selection 



Filter 


Channel 


Central A (/im) 


Width (/im) 


[Fe II] 


Blue 


1.644 


0.0247 


CH4 H 1% L 


Red 


1.628 


0.0174 


H2 1-0 8(1) 


Blue 


2.124 


0.0261 


Br-gamma 


Blue 


2.169 


0.0295 


Kcont 


Red 


2.272 


0.0352 
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Table 2. Velocities and Ejection Epochs 



Region Projected Vel Origin 
Ap. Number (km s"'^) Epoch'' 



Doppler Vel Total Vel 0^ 
(km s-i) (km s"!) (deg) 



Distance'' Origin i^" 
(AU) Epoch'= (deg) 



A; Western Arc 
1 
2 

B: SW Region 
3 
4 

5 



9 

10 
11 
12 

C: S Clump 
13 
14 
15 

16 
17 
18 

D: Southern Arc 
19 
20 

21 
22 
23 
24 
25 



87.1±3.4 
57.2±5.8 

81.2±5.2 

145.7±3.7 
145.4±2.2 
159.8±2.6 
164.2±2.5 
144.2±2.0 
138.2±1.9 
122.7±1.1 
155.4±1.3 
184.6±2.2 

101.1±2.0 
100.2±1.8 
91.0±3.5 

68.9±2.6 
96.5±1.4 
119.4±2.0 

116.7±2.0 
135.7±1.7 

150.3±1.2 
103.2±0.7 
98.2±0.8 
86.4±1.6 
77.5±1.9 



E: SE Clump and Filament 
26 61.1±1.1 



1902±8 
1830±36 

1848±20 
1906±5 
1911±3 
1917±3 
1903±3 
1877±4 
1865±4 
1852±3 
1871±2 
1883±3 

1854±6 
1885±4 
1902±8 

1901±8 
1898±3 
1897±4 

1898±4 
1907±3 
1910±2 
1887±2 
1895±2 
1887±4 
1900±5 

1901±4 



-52±5 



96.7±6.0 



38±5 



107.3±5.3 



38±5 

38±5 
38±5 
38±5 



125.3±5.4 

122.7±5.4 
140.8±5.3 
155.0±5.1 



-32.9±3.4 



3286 



1845±17 



21.0±3.7 



2807 



1883±10 



17.7±3.5 

18.0±3.5 
15.6±3.2 
14.9±3.0 



2945 

2877 
3032 
3236 



1896±8 

1896±8 
1905±7 
1908±6 



-13 
-13 
-17 
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Table 2 — Continued 



Region 


Projected Vel 


Origin 


Doppler Vel 


Total Vel 


0b 


Distance'^ 


Origin 


c 


Ap. Number 


(km s-i) 


Epoch^- 


(km s-l) 


(km s-i) 


(deg) 


(AU) 


Epoch-^ 


(deg) 


27 


44.3±1.7 


1892±9 














28 


29.2±1.6 


1878±14 


48±5 


56.2±5.3 


58.7±2.5 


1549 


1876±14 


+18 


29 


33.6±2.2 


1881±16 


48±5 


58.6±5.4 


55.1±3.0 


1608 


1877±15 


+17 


30 


39.9±1.9 


1874±12 


48±5 


62.4±5.3 


50.0±2.9 


1796 


1870±14 


+8 


31 


43.8±2.2 


1879±12 


48±5 


65.0±5.4 


47.6±3.1 


1773 


1878±13 


+3 


32 


52.6±1.4 


1886±7 


48±5 


71.2±5.2 


42.4±2.9 


1837 


1885±11 


-1 


33 


52.3±1.9 


1862±11 


48±5 


71.0±5.4 


42.5±3.1 


2214 


1859±14 


-1 


F; SE Arc 


















34 


57.5±1.4 


1869±7 


48±5 


74.9±5.2 


39.9±2.9 


2218 


1866±13 


-12 


35 


63.3±1.3 


1848±7 


48±5 


79.4±5.2 


37.2±3.0 


2685 


1846±14 


-1 


36 


70.4±1.7 


1852±7 


48±5 


85.2±5.3 


34.3±3.2 


2811 


1850±13 


-7 


37 


80.7±1.5 


1857±6 


48±5 


93.9±5.2 


30.7±3.1 


3022 


1854±12 


-13 


38 


80.2±2.9 


1825±13 


48±5 


93.5±5.8 


30.9±3.6 


3618 


1822±16 


—11 


39 


82.9±3.0 


1817±14 


48±5 


95.8±5.8 


30.1±3.6 


3880 


1814±17 


-5 


G: East Edge Connecting SE 


Arc and NE Region 












40 


95.0±1.1 


1838±4 














H: NE Clump 


















41 


116.1±1.1 


1863±3 














42 


78.8±1.4 


1830±6 


58±5 


97.8±5.2 


36.4±2.5 


3685 


1827±13 


+10 


43 


85.3±2.5 


1862±8 


58±5 


103.2±5.6 


34.2±2.9 


3199 


1860±11 


+10 


44 


60.5±2.6 


1817±16 


58±5 


83.8±5.6 


43.8±2.9 


3377 


1815±16 


+14 


45 


58.5±1.1 


1839±6 


58±5 


82.4±5.1 


44.8±2.2 


2949 


1837±13 


+19 


46 


58.4±0.9 


1843±5 


58±5 


82.3±5.1 


44.8±2.1 


2884 


1840±13 


+22 


47 


61.9±4.3 


1867±19 














48 


42.1±6.0 


1827±51 














I: South side of NE Region 
















49 


86.2±1.7 


1848±6 














50 


74.7±1.1 


1832±5 


58±5 


94.5±5.1 


37.8±2.4 


3521 


1829±13 


+23 


51 


74.5±1.3 


1833±6 


58±5 


94.4±5.2 


37.9±2.4 


3497 


1830±13 


+19 


52 


82.9±1.7 


1829±7 














53 


77.5±1.2 


1802±6 
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Table 2 — Continued 



Region 


Projected Vel 


Origin 


Doppler Vel 


Total Vel 




Distance'' 


Origin 


C 


Ap. Number 


(km s-i) 


Epoch" 


(km s-l) 


(km s-i) 


(deg) 


(AU) 


Epoch"= 


(deg) 


54 


120.9±1.8 
















55 


125.8±2.0 


1845±5 














56 


128.8±1.3 


1840±3 














57 


153.2±3.0 


1877±5 














J: North side of NE Region 
















58 


110.3±2.9 


1844±8 














59 


128.3±2.3 


1852±6 














60 


107.3±2.5 


1843±7 














61 


94.4±2.2 


1835±8 














62 


83.3±2.9 


1831±12 














63 


111.7±1.2 


1861±3 














64 


91.3±0.7 


1843±2 














65 


85.3±1.6 


1850±6 


-62±5 


105.5±5.3 


-36.0±2.4 


3524 


1848±14 


-14 


66 


87.9±2.0 


1862±6 


-62±5 


107.6±5.4 


-35.1±2.5 


3316 


1860±13 


-12 


67 


96.7±5.0 


1869±14 


-62±5 


114.9±7.1 


-32.7±3.6 


3388 


1867±16 


-8 



^From the total projected motion in the plane of the sky. 
''Orientation angle with respect to the plane of the sky. 
°Prom the total space velocity. 

''Distance from the star in AU at the inclination angle. 
"The angle or latitude relative to the equatorial plane. 



Table A. Positions and Proper Motions 



Aperture Dist. from PA'' Relative Proper Motion Avg Homunculus'^ True Motion*^ 

Number Star (arcsec) (deg) RA (mas/yr) Dec (mas/yr) N** Sig'^ Exp. (mas/yr) RA (mas/yr) Dec (mas/yr) 



1 


0.86 


-77.29 


-2.56±0.58 


0.46±0.17 


8 


0.93 


5.39 


-7.82±0.58 


1.64±0.17 


2 


0.94 


-80.21 


0.60±0.78 


-0.63±0.69 


6 


0.93 


5.92 


-5.23±0.78 


0.38±0.69 


3 


1.20 


-99.73 


0.77±0.12 


-2.09±0.92 


8 


0.84 


7.53 


-6.65±0.12 


-3.36±0.92 


4 


1.38 


-124.15 


-2.78±0.35 


-4.01±0.56 


7 


0.90 


8.70 


-9.98±0.35 


-8.89±0.56 


5 


1.30 


-127.36 


-2.45±0.13 


-4.90±0.38 


7 


0.92 


8.19 


-8.96±0.13 


-9.88±0.38 


6 


1.35 


-130.02 


-2.98±0.18 


-5.71±0.43 


7 


0.91 


8.50 


-9.49±0.18 


-11.17±0.43 


7 


1.59 


-142.85 


-3.43±0.37 


-3.71±0.25 


7 


0.89 


10.03 


-9.48±0.37 


-11.71±0.25 


8 


1.75 


-147.30 


0.05±0.15 


-2.58±0.33 


8 


0.92 


11.00 


-5.90±0.15 


-11.84±0.33 


9 


1.83 


-149.59 


0.40±0.20 


-1.53±0.27 


7 


0.91 


11.52 


-5.43±0.20 


-11.46±0.27 


10 


1.77 


-153.81 


0.15±0.12 


-0.20±0.16 


7 


0.89 


11.15 


-4.77±0.12 


-10.20±0.16 


11 


1.97 


-157.71 


-0.34±0.10 


-1.90±0.20 


9 


0.91 


12.37 


-5.03±0.10 


-13.34±0.20 


12 


2.14 


-161.66 


-3.36±0.38 


-2.36±0.07 


7 


0.85 


13.46 


-7.59±0.38 


-15.14±0.07 


13 


1.44 


-165.25 


0.42±0.17 


-0.34±0.33 


7 


0.90 


9.04 


-1.88±0.17 


-9.08±0.33 


14 


1.14 


-166.95 


0.13±0.05 


-2.09±0.31 


7 


0.93 


7.16 


-1.49±0.05 


-9.07±0.31 


15 


0.89 


-172.76 


1.21±0.10 


-2.77±0.62 


7 


0.94 


5.60 


0.51±0.10 


-8.33±0.62 


16 


0.68 


-176.65 


-0.53±0.41 


-1.98±0.23 


6 


0.61 


4.30 


-0.78±0.41 


-6.27±0.23 


17 


0.99 


-172.40 


0.78±0.16 


-2.71±0.19 


8 


0.93 


6.20 


-0.04±0.16 


-8.85±0.19 


18 


1.22 


-169.59 


-0.36±0.14 


-3.24±0.33 


7 


0.94 


7.69 


-1.75±0.14 


-10.81±0.33 


19 


1.19 


-178.95 


-0.66±0.17 


-3.20±0.33 


8 


0.92 


7.48 


-0.79±0.17 


-10.68±0.33 


20 


1.27 


-175.76 


-1.16±0.10 


-4.33±0.29 


8 


0.93 


8.02 


-1.75±0.10 


-12.33±0.29 


21 


1.36 


178.64 


-1.24±0.08 


-5.17±0.20 


7 


0.92 


8.58 


-1.03±0.08 


-13.75±0.20 


22 


1.16 


172.11 


0.24±0.07 


-2.17±0.09 


6 


0.93 


7.30 


1.24±0.07 


-9.39±0.09 


23 


1.03 


172.10 


0.56±0.04 


-2.50±0.14 


6 


0.92 


6.45 


1.44±0.04 


-8.89±0.14 


24 


0.97 


165.09 


0.58±0.13 


-1.73±0.26 


7 


0.89 


6.10 


2.15±0.13 


-7.63±0.26 


25 


0.78 


164.06 


0.34±0.24 


-2.20±0.23 


9 


0.88 


4.89 


1.68±0.24 


-6.91±0.23 


26 


0.61 


161.24 


-0.83±0.12 


-1.97±0.15 


7 


0.92 


3.82 


0.40±0.12 


-5.59±0.15 


27 


0.48 


158.17 


-0.59±0.08 


-1.25±0.29 


7 


0.96 


3.00 


0.53±0.08 


-4.03±0.29 


28 


0.35 


156.25 


-0.97±0.24 


-0.67±0.16 


7 


0.91 


2.20 


-0.08±0.24 


-2.68±0.16 


29 


0.40 


140.81 


-1.64±0.37 


-1.15±0.11 


7 


0.87 


2.49 


-0.07±0.37 


-3.08±0.11 


30 


0.50 


136.63 


-0.95±0.26 


-1.19±0.22 


7 


0.91 


3.12 


1.19±0.26 


-3.46±0.22 


31 


0.52 


120.74 


0.80±0.34 


-0.06±0.18 


6 


0.87 


3.28 


3.62±0.34 


-1.74±0.18 
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Table A — Continued 



)erture 


Dist. from 




Relative Proper Motion 




Avg 


Homunculus'' 


True Motion'^ 


umber 


Star (arcsec) 


(deg) 


RA (mas/yr) 


Dec (mas/yr) 




Sig« 


Exp. (mas/yr) 


RA (mas/yr) 


Dec (mas/yr) 


32 


0.59 


118.35 


1.28±0.19 


0.20±0.17 


7 


0.90 


3.73 


4.56±0.19 


-1.57±0.17 


33 


0.71 


113.46 


0.38±0.31 


-0.01±0.15 


7 


0.86 


4.44 


4.46±0.31 


-1.78±0.15 


34 


0.74 


135.01 


1.30±0.23 


-1.61±0.10 


8 


0.95 


4.64 


1.98±0.23 


-4.89±0.10 


35 


0.93 


131.20 


-3.24±0.22 


-1.82±0.10 


7 


0.96 


5.87 


1.18±0.22 


-5.69±0.10 


36 


1.01 


130.18 


-2.65±0.24 


-1.96±0.19 


7 


0.94 


6.38 


2.22±0.24 


-6.07±0.19 


37 


1.13 


129.86 


-1.59±0.24 


-1.77±0.15 


7 


0.90 


7.09 


3.85±0.24 


-6.32±0.15 


38 


1.35 


129.26 


-3.49±0.51 


-1.30±0.10 


7 


0.89 


8.50 


3.09d=0.51 


-6.68±0.10 


39 


1.46 


126.55 


-6.47±0.51 


-2.08±0.16 


8 


0.86 


9.18 


0.91±0.51 


-7.55±0.16 


40 


1.49 


100.34 


-0.51±0.08 


1.01±0.17 


6 


0.84 


9.35 


8.69±0.08 


-0.67±0.17 


41 


1.55 


70.35 


0.70±0.20 


0.66±0.06 


7 


0.87 


9.77 


9.90±0.20 


3.94±0.06 


42 


1.29 


66.91 


-0.72±0.23 


-0.66±0.09 


7 


0.92 


8.14 


6.77±0.23 


2.53±0.09 


43 


1.15 


65.23 


0.92±0.42 


-0.79±0.16 


6 


0.93 


7.25 


7.50±0.42 


2.24±0.16 


44 


1.06 


71.62 


-0.82±0.46 


-1.58±0.09 


7 


0.80 


6.69 


5.53±0.46 


0.53±0.09 


45 


0.91 


62.92 


0.07±0.14 


-1.17±0.13 


7 


0.96 


5.73 


5.17±0.14 


1.44±0.13 


46 


0.89 


56.25 


0.19±0.14 


-0.80±0.06 


7 


0.95 


5.59 


4.84±0.14 


2.30±0.06 


47 


0.81 


60.32 


-0.21±0.10 


1.29±0.76 


7 


0.83 


5.08 


4.21±0.10 


3.81±0.76 


48 


0.70 


45.56 


0.17±0.41 


-1.14±1.00 


7 


0.96 


4.43 


3.33±0.41 


1.96±1.00 


49 


1.27 


42.79 


1.71±0.19 


-2.49±0.24 


7 


0.86 


8.01 


7.15±0.19 


3.39±0.24 


50 


1.21 


46.76 


0.38±0.19 


-1.82±0.08 


8 


0.86 


7.63 


5.94±0.19 


3.41±0.08 


51 


1.20 


53.96 


-0.02±0.22 


-1.36±0.07 


7 


0.87 


7.58 


6.10±0.22 


3.10±0.07 


52 


1.37 


57.44 


-0.70±0.20 


-0.81±0.24 


8 


0.94 


8.62 


6.57±0.20 


3.83±0.24 


53 


1.48 


53.85 


-2.23±0.18 


-0.69±0.12 


8 


0.89 


9.28 


5.26±0.18 


4.79±0.12 


54 


1.66 


52.73 


-0.21±0.23 


1.29±0.23 


7 


0.83 


10.41 


8.08±0.23 


7.60±0.23 


55 


1.89 


59.35 


-0.74±0.28 


0.51±0.21 


7 


0.83 


11.88 


9.49±0.28 


6.57±0.21 


56 


2.00 


60.32 


-1.02±0.19 


0.25±0.15 


8 


0.85 


12.55 


9.89±0.19 


6.47±0.15 


57 


1.86 


51.14 


2.70±0.48 


0.33±0.21 


6 


0.85 


11.67 


11.79±0.48 


7.66±0.22 


58 


1.67 


38.10 


0.47±0.16 


-0.88±0.49 


7 


0.85 


10.48 


6.94±0.16 


7.37±0.49 


59 


1.85 


37.22 


0.89±0.30 


-0.56±0.29 


7 


0.84 


11.63 


7.93±0.30 


8.70±0.29 


60 


1.65 


33.17 


-0.68±0.24 


-0.12±0.37 


7 


0.88 


10.36 


4.99±0.24 


8.55±0.37 


61 


1.51 


32.52 


-0.94±0.23 


-0.40±0.33 


7 


0.92 


9.49 


4.16±0.23 


7.60±0.33 


62 


1.36 


32.07 


-0.33±0.33 


-0.85±0.42 


8 


0.93 


8.53 


4.21±0.33 


6.38±0.42 
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Table A — Continued 



Aperture 


Dist. from 




Relative Proper Motion 




Avg 


Homunculus'^ 


True Motion'' 


Niiiiilx^r 






HA (inas/yr) 


Dec (luas/yr) 


N'' 


Sis" 


Ex]>, (uias/\'r) 


HA. (iiias/\a.) 


Dec (iiias/\a") 


63 


1.52 


24.55 


-1.40±0.16 


1.24±0.15 


8 


0.91 


9.54 


2.56±0.16 


9.92±0.15 


64 


1.39 


26.87 


-1.63±0.08 


0.23±0.10 


7 


0.92 


8.76 


2.33±0.08 


8.05±0.10 


65 


1.24 


24.87 


-0.74±0.27 


0.32±0.13 


7 


0.94 


7.80 


2.54±0.27 


7.40±0.13 


66 


1.18 


22.27 


0.05±0.28 


0.64±0.22 


8 


0.93 


7.45 


2.88±0.28 


7.53±0.22 


67 


1.24 


18.51 


0.18±0.31 


1.05±0.84 


7 


0.92 


7.82 


2.66±0.31 


8.46±0.84 



^Position angle measured from north through east. 
''Number of cells included in the aperture. 

'^The weighted average of the significance of the cells included in the aperture. 
■^The radial expansion of the Homunculus at this position from the star. 
^The relative motion corrected for the expansion of the Homunculus. 
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Fig. 1. — The continuum (top row) and continuum-subtracted (bottom row) H2, Br7, and 
[Fe II] images. The grey scale is arbitrary. All of the images are oriented north at the top 
and east to the left. 




Fig. 2a. — The combined Br7, H2, [Fe II] image without the continuum subtraction; Br7 
(red), H2 (green), and [Fe II] (blue). North is up and east is left. 
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Fig. 2b. — The combined line - continuum image; Br7 (red), H2 (green), and [Fe II] (blue). 
This figure shows the inner part of the Homunculus and region of the butterfiy. Note the very 
different distribution of the [Fe II] emission compared with Br7 and H2. Br7 and H2 have 
nearly identical distributions, therefore their combined image appears white or yellowish. 
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Fig. 3.— The NICMOS 2.15/im (ISmith fc Gehrz 



20001 1 (top left) compared with the NICI 



i^-band continuum i mage (top right) and on the bottom row the central region from the 



NaCO broadband (IChesneau et al. 



1200,^ (left) and the NICI K-haxid images (right). All 
four images have been matched to the same spatial scale and orientation, north at the top 
and east to the left. The field of view, 16.5", is the same. A different contrast and intensity 



scaling factor has been used for the upper and lower NICI images 
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Fig. 4. — The normalized fluxes vs distance from the star for knot C' (left panel) and knot 
D' (right panel) measured at different times: in order, from the upper solid curve, 2002.88, 
2003.95, 2004.99 to 2005.98, the lower dash-dot-dash line. 
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Fig. 5. — The scaling factor between the NaCO and NICI images as a function of time. Using 
the measured scaling factors with the known pixel scales, we find that the Homunculus lobes 
were ejected in 1850±4.0. 
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2 1 0-1-2 

field position (") 

Fig. 6. — Plot of the relative, co-moving vector motions between the NICI (2009.1) and the 
NaCO (2002.9) Ks-continuum images. The length of the arrows is proportional to the size 
of the total motion. The arrows are color coded by the significance of the correlation (see 
text); red, 0.8 to 1.0; blue, 0.6 to 0.8 and green, 0.4 to 0.6. The circle in the upper left is 
the size of the correlation region. Also see Figure 8. 



Fig. 7. — Example of a filamentary structure showing a significant motion between NaCO 
Ks (left) and NICI (middle) images. The superposition of the NICI and NaCO contour plots 
(right) illustrates the motion of the filamentary structure. The South-Eastern tip (at x=1.2, 
y =-1.0) of the filament has a tranverse motion of about 8 mas/yr while the bright knot (at 
x=0.8, y =-0.7) has a smaller (about 4 mas/yr) tranverse motion in the same direction. 
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Fig. 8. — A map of the apertures and the vectors for the true or intrinsic motions in the 
fixed x,y grid (see Appendix A) plotted on the NICI K-band continuum image. The arrow 
lengths are proportional to the total motion. 
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1820 1840 1860 1880 1900 1920 
Ejection Date 

Fig. 9. — The distribution of the ejection epochs from the motions in Table 2. The brackets 
indicate the times of the great eruption (1838 - 1858) and the second eruption (1887 - 1895). 
The mean uncertainty in the ejection epochs is 7 years. Most of the apertures have an 
uncertainty of 5 to 10 years. The errors are quoted in Table 2. 
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Fig. 10. — A map of the larger regions and apertures (see Table 2) plotted on the NICI 
K-band continuum image. The parallel dashed lines from lower left to upper right note the 
locations of the slits used to measure radial velocities from Smith (2006). 
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Fig. 11. — The inclination i angle (40.7°) between our line of sight and the axis of the 
bipolar Homunculus, adapted from Davidson et al. (2001). This agrees with the conventional 
definition of i for binary star orbits, if we substitute the Homunculus equatorial plane for the 
orbit plane. The total polar diameter of the Homunculus is roughly 40000 AU, much larger 
than the area shown in Fig. 11. Caveat: In fact the axial symmetry is only an approximation. 
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Fig. 12. — "Side view" of the aperture positions listed in Table 2 and shown in Fig. 8. These 
are 3-D locations projected into the xz plane defined in Fig. 11, but the area shown here is 
much smaller. Since this is a projected view rather than a cross-section, random samples of 
a circular or toroidal configuration would create a band connecting the upper left and lower 
right sets of points; so the distribution plotted here does not provide evidence for such a 
structure. 



